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ABSTRACT 

Context: It has been recently shown that GRB 090618, observed by AGILE, Coronas Photon, Fermi, Konus, Suzaku and Swift, is 
composed of two very different components: episode 1, lasting 50 s, shows a thermal plus power-law spectrum with a characteristic 
temperature evolving in time as a power law; episode 2 (the remaining 100 s) is a canonical long GRB. We have associated episode 1 
to the progenitor of a collapsing bare core leading to the formation of a black hole: what was defined as a "proto black hole". 
Aims: In precise analogy with GRB 090618 we aim to analyze the 89s of the emission of GRB 101023, observed by Fermi, Gemini, 
Konus and Swift, to see if there are two different episodes: the first one presenting a characteristic black-body temperature evolving 
in time as a broken power law, and the second one consistent with a canonical GRB. 

Methods: To obtain information on the spectra, we analyzed the data provided by the GBM detector onboard the Fermi satellite, and 
we used the heasoft package XSPEC and RMFIT to obtain their spectral distribution. We also used the numerical code GRBsim to 
simulate the emission in the context of the fireshell scenario for episode 2. 

Results: We confirm that the first episode can be well fit by a black body plus power-law spectral model. The temperature changes with 
time following a broken power law, and the photon index of the power-law component presents a soft-to-hard evolution. We estimate 
that the radius of this source increases with time with a velocity of 1.5 x I0*km/s. The second episode appears to be a canonical GRB. 
By using the Amati and the Atteia relations, we determined the cosmological redshift, z ~ 0.9 ± 0.0S4(stat.) ± 0.2(sys.). The results 
of GRB 090618 are compared and contrasted with the results of GRB 101023. Particularly striking is the scaling law of the soft X-ray 
component of the afterglow. 

Conclusions: We identify GRB 090618 and GRB 101023 with a new family of GRBs related to a single core collapse and presenting 
two astrophysical components: a first one related to the proto-black hole prior to the process of gravitational collapse (episode 1), 
and a second one, which is the canonical GRB (episode 2) emitted during the formation of the black hole. For the first time we are 
witnessing the process of a black hole formation from the instants preceding the gravitational collapse up to the GRB emission. This 
analysis indicates progress towards developing a GRB distance indicator based on understanding the P-GRB and the prompt emission, 
as well as the soft X-ray behavior of the late afterglow. 

Key words. Gamma-ray burst: individual: GRB 101023 — Black hole physics 



1. Introduction 



Discovered at the end of the 60s ("Strong 1975[ l, gamma-ray 
bursts (GRBs) are extremely intense flashes of hard X-radiation, 
coming from random directions in the sky at unpredictable times 
and typically lasting from a fraction of a second up to a few min- 
utes. They are detected by satellites in low Earth orbit at a rate 
of ~0.8 events/day. As outlined by breakthrough observations in 
the last ~15 years, these phenomena are by far the most ener- 
getic sources in the Universe, observed in a range of cosmologi- 
cal redshift 0.0084 < z < 9 (|Salvaterra et al.||2^[Tinvir et al.| 



2009 Cucchiara et al. 201 It , with isotropic equivalent radiated 
energy Eiso in the range lO'* - 10^^ erg and a th eoretically pre - 



dicted upper limit to their energies of 10'' erg ([Ruffini 



2011 



Since the early observation by BATSE (Meegan et al. 1992 1, 
they have been divided into two classes: the short GRBs, with a 
characteristic duration of Tgo < 2 s, and the long GRBs, with a 



characteristic Tgo > 2 s ( 


Dezalayetal. 1992| Klebesadel 1992 


IKouvehotou et al.||1993 





Analysis of the GRBs within the fireshell model (see e.g. 
IRuffini et al.| |200H [2009 1 and references therein) has led to 
identifying a canonical GRB structure described by two parame- 
ters: the total energy E^^^ of the initially optically thick electron- 
positron plasma and its baryon load B = Mbc^ IE"gj. To this 
information characterizing the source is added the information 
on the density and filamentary distribution of the circumburst 
medium (C BM) ( .Ruflini et al.[ |2004b, .20"05) |Patricelli et al] 
i20T0l|20TT] l. 

Within this model the structure of a canonical GRB has been 
identified. It is composed by a proper-GRB (P-GRB), followed 
by an extended afterglow. The P-GRB originates at the mo- 
ment of transparency of the relativistically expanding electron- 
positron plasma. The extended afterglow originates in the colli- 
sion of the ultra-relativistic baryons with the filamentary struc- 
ture of the CBM. The acceleration process of the baryons oc- 
curs in the optically thick phase of the self-accelerating electron- 
positron plasma. This explains the spiky emission observed in 
the prompt radiation ( |RufRni et al.| |2002[ ). The average density. 
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the porosity, and the dimensions of the clouds in the CBM are in 
turn determined (see e.g. 'Ruffini et al. 2006 Bernardini et al. 



2007 ; Dainotti et al.| [200 7 ; Ca ito et al.||2009f . 

This model has allowed the nature of long GRBs to be ex 
plained and two new classes of short bursts to be introduced. A 
first class contains the disguis ed short GRBs ([B ernardini et al 



2007) [Caito et al.j [2009] pmO\ [De Barros etaL||2011 ): just long 
GRBs exploding in low density CBM (n = 10'' part/ cirP), and 



often referred to as short GRBs in the literature (see e.g. Gehrels 



et al. 2005 i. A second class contains the genuine short GRBs, 



theoretically foreseen in [Ruffini et al. (2001 ) as canonical GRBs 
occurring in the limit of a very low baryon load, B < 10'^. This 
new class of genuine short GRBs is expected to occur on a much 
shorter time scale, Tgo < 10"^ - 10"^ s. 

With the observation of GRB 090618, a novel situation has 
occurred with respect to the above classification. It had been 
shown in the pioneering works of Felix Ryde and his collabora- 
tors (Ryde 2004) that, in the early emission of selected BATSE 
sources and also in some Fermi sources, a characteristic thermal 
component is present with temperature changing in time follow- 
ing a broken power law ( Ryde 2004, 2005, Ryde & Pe'er , 2009 ). 
They attempted to interpret this emission within the GRB fireball 
model (see e.g. Pe'er et al. 2007 1. 



Ruffini et al. ( 2010a i showed that two very different episodes 



occur in GRB 090618: episodes 1 and 2. Episode 1 presents an 
emission "a la Ryde". There it was proposed that such an emis- 
sion, alternatively to the Ryde interpretation, had to be inter- 
preted as originating in a new kind of source in the late phase of 
a core collapse. The concept of proto-black hole was introduced 
there. Episode 2 was shown to be consistent with a canonical 
GRB. 

Details of the data analysis showing the characteristic broken 
power law temporal variation of the temperature of the thermal 
component of episode 1 are presented in|Izzo et al.|((201 1 ). The 
radius of the emitting region and its time variation have been 
determined as well, along with the details of the GRB emis- 
sion of episode 2, including the P-GRB structure, the porosity 
of the interstellar medium, the baryon load B, and the total en- 
ergy. Identifying these two components has been made possible 
by the extraordinary coincidence of three major factors for this 
GRB: 1) precise determination of the cosmological redshift of 
this source z = 0.54, implying the fortunate occurrence of a very 
close source with an energy Ei,,, - 2.7 x lO^^* erg; 2) joint ob- 
servations by several X and gamma-ray telescopes; 3) the ex- 
ceptional dataset on the instantaneous spectral distribution, light 
curve, and luminosity variation of this source (see section[2]). 

There is a striking morphological analogy between GRB 
101023 and GRB 090618 (see Figs.|2]and[T]l. Both light curves 
present a first emission that lasts ~ 50 s, fo lowed by a spikier 
structure in the remaining part. We identify the first 45 s of GRB 
101023 with episode 1 and the remaining 44 s with episode 2 
(a canonical GRB). There is, however, a substantial difference 
between these two sources. In the present source, GRB 101023, 
the cosmological redshift is unknown. This has not been a draw- 
back for us but a challenge that probes our understanding of the 
GRB phenomenon. It is interesting, as a rough estimate, that if 
one were to assume that the two sources, GRB 101023 and GRB 
090618 had not only the same morphology but also the same en- 
ergy Eiso, one would infer z = 1 for the cosmological redshift 
of GRB 101023. A main result of this article is that, assuming 
the validity of the Amati relation (see 



Amati et al. 



P009| and 



references therein) and Atteia criteria (Atteia 20031, it is pos- 
sible to theoretically derive an expected cosmological redshift 
z = 0.9 + 0.084(5taf.) ± Q.2(sys.) for episode 2. 



What is most striking is that we can have an independent ver- 
ification of this redshift by comparing the late part of the after- 
glows of the two sources. Since we have verified that both GRB 
090618 and GRB 101023 have similar energetics, and under the 
hypothesis of the same progenitor mechanism, we compare and 
contrast the luminosities of both GRBs in the late X-ray after- 
glow emission. We know that the X-ray afterglow is related to 
the residual kinetic energy of the outflow, although we do not 
attempt here to present a theoretical model for this emission. We 
rescaled, in the observed time interval and energy range, the X- 
ray afterglow luminosity of GRB 090618 for different redshifts 



in an interval between 0.04 < z < 3 (see Fig. 15 1. The striking 
coincidence for z = 0.9 is presented in Fig. [T4| 

In section [2] we summarize the results of GRB 090618 and 
identify episode 1 and episode 2. In section [3] we present the 
observations of GRB 101023 by the different satellites. In sec- 
tion [4[ we give a brief summary of the fireshell scenario. In sec- 
tion p] we perform a spectral analysis of episodes 1 and 2 of this 
GRB. In section [6] we try to identify the P-GRB of the gamma- 
ray burst, taking different time intervals into account along the 
entire emission. In section [7] we present the methods we used 
to constrain the redshift. In section [8] after interpreting the sec- 
ond episode as a canonical GRB within the fireshell model, we 
build its light curve and spectrum. In Section [9| we go into fur- 
ther detail in the analysis of the first episode, making clear the 
evolution of the thermal component and the radius of the out- 
ermost shell and establishing the complete correspondence with 
GRB 090618. Finally, in Section 10 we present the conclusions. 



2. Brief summary of GRB090618 analysis 

We recall that GRB 0908618 is one of the most energetic among 
the nearest sources, with an isotropic energy of Ejso - 2.7 x 10^"* 
erg, at redshift z= 0.54. It has been observed in a wide energy 
range by many satellites, such as as Fermi GBM (Meegan et al. 



2009), Swift-BAT ( [Gehrels et ai:][2009i), AGILE ( ,Longo et al. 
200 9|, Konus-WIN D ( [Golenteskii et al.] [2009l l, Suzaku-WAJV 



(Ko no et al.j [2009] ), and CORONAS -PHOTON (Kotov et al. 



20081, and by many onground telescopes. We have shown (see 
the work of [Izzo et al.[ [201 l| l that the light curve is quite par- 
ticular, as it consists of two different emissions, of 50 s and 100 
s of duration. A time-resolved spectral analysis showed that the 
first part is well fit by a black body and an extra power-law com- 
ponent. The temperature decays with time following a broken 
power law, in agreement with the results found by Ryde and 
collaborators ( [Ryde & Pe'er[ [2009| l. The first power law has 
an index a^T - -0.33 + 0.07, and the second one has an in- 
dex bkT = -0.57 ± 0.11. The evolution of the radius of the 
black body emitter has also been studied, finding an initial ra- 
dius of 12000 km, expanding in the early phase with a veloc- 
ity of ~ 4000 km/s. By analyzing it within the fireshell model, 
we concluded that the first episode cannot be either a GRB or 
part of a GRB. Indeed, we relate this episode to the phases just 
preceding the gravitational collapse and define it as a "proto- 
black hole": the latest phase of the collapsing bare core leading 
to the black hole formation and the simultaneous emission of 
the GRB (Ruffini et al. 2010a| . In this interpretation, the radius 
r„„ only depends on the observed energy flux of the black body 
component (pobs^ the temperature kT and the luminosity distance 
to the source D. Episode 2 was identified as a canonical GRB, 
which comes from the black hole formation process. The first 
4 s were identified as the P-GRB, and its spectrum is well fit 
by a black body with an extra power-law component, the lat- 
ter mainly caused by the early emission of the extended after- 
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Fig. 1. Count light curve of GRB 090618 obtained from Fermi 
GBM detector, with a bin time of 1 s, and showing two-episode 
nature of the GRB. 

glow. We found a P-GRB temperature of kT = 29.22 + 2.21 
keV and a dyadosphere energy of the whole second episode 
of E^g, = 2.49 X lO^-' erg. We performed a numerical simula- 
tion with the numerical code GRBsim and found a baryon load 
B - (1.98±0.15)x lO^-'and a Lorentz Gamma factor at the trans- 
parency of r = 495 + 40. From this analysis we concluded that 
we are in the presence of a very interesting source, because for 
the first time we can witness the process of formation of a black 
hole from the phases just preceding the gravitational collapse to 
the GRB emission. 



3. Observations of GRB 101023 



On 23 October 2010 the Fermi GBM ( |GCN circular) [20T0l l de- 
tector was triggered by a source quite similar to GRB 090618, 
with a trigger time of 309567006.72696 8 (in MET seco nds) 
The burst was also detected by BAT (Saxton et al. 2010| l (i 
Fig. [3]l, onboard the Swift satellite (Gehrels et al. ' 



2004K 



see 
with 

a trigger time of 436981 (in MET seconds) and the following 
location coordinates: RA{J2000) = 2lhlhn49s, DeciJlQQQ) = 
-65°23'37" with an uncertainty of 3 arcmin. The Swift-XRT 
detector ( jPage et aL 2010[ Burrows et al. 2005 i has also ob- 
served this source from 88 s to 6.0 ks after the BAT trigger. 
GRB 101023 was also detected by the Wind instrument onboard 
Konus satellite, in the energy range (10 - 770) keV (Gol enetskii| 
|et al.|[201 0). The inferred location is in complete agreement with 
that determined by Swift and Fermi. Moreover, there have been 
detections in the optical band by the Gemini telescope (Levan et 
[aTl lSoTOl l. 

The GBM light curve (Fig. |2]i shows two major pulses. The 
first one starts at the trigger time and lasts 45 s. It consists of 
a small peak that lasts about 10 s, followed by a higher emis- 
sion that decays slowly with time. The duration, as well as the 
topology of this curve, lead us to think that this may not be a 
canonical GRB, but its origin may lie on another kind of source, 
which remains unidentified. The second pulse starts at 45 s af- 
ter the trigger time and lasts 44 s. It presents a peaky structure, 
composed of a short and weak peak at the beginning, followed 
by several bumps, big not only in magnitude but also in duration. 
This second emission, in contrast, does have all the characteris- 
tics that describe a canonical GRB ( jRuffini et al.||20I0cl l. 



qiq-cif)M:-n2-lifiiaifl2395l->iQQ4ifi[i: glg-i:ape[:-n2-bn1Q1Q33951-WM-lu 
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Fig. 2. Count Hght curve of GRB 101023 obtained from the 
Fermi GBM detector, with a bin time of 1 s. The time is given 
with respect to the GBM trigger time of 22:50:04.73 UT, 2010 
October 23. The plot was obtained with the RMFIT program. 
The two-episode nature of the GRB is shown in analogy with 
GRB 090618. 
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Fig. 3. Count light curve of GRB 101023 obtained from the 
Swift XRT detector. 



4. Theoretical model considered: Fireshell Scenario 

In the fireshell scenario, the GRB emission comes from a pro- 
cess of vacuum polarization, resulting in pair creation in the 
so-called dyadosphere. In the process of gravitational collapse 
to a black hole (Ruffi ni et al.||201 0b), an e"^ plasma is formed 
in thermal equilibrium, with total energy E'j^,. The annihilation 
of these e"^ pairs occurs gradually and is confined in a shell, 
called "fireshell". This shell self-accelerates to relativistic ve- 
locities, engulfing the baryonic matter (of mass Mb) left over in 
the process of collapse and reaching a thermal equilibrium with 
it (Ruflini et al. 20001. The baryon loading is measured by the 
dimensionless parameter B = Msc^/E^g,. The fireshell continues 



to self-accelerate up to relativistic velocities ( |Ruflini|[l999| l until 
it reaches the transparency condition. At this time we have a first 
flash of radiation, the P-GRB ( [Ruffini et aL|[200T) . The energy 
released in the P-GRB is a fraction of the initial energy of the 



3 



A. V. Penacchioni et al.: A double astrophysical component in GRB 101023 



dyadosphere isf^,. The residual plasma of leptons and baryons 
interacts with the circumburst medium (CBM) as it expands, 
giving rise to multi-wavelength emission: the "extended" after- 
glow. However, owing to these collisions, the plasma starts to 
slow down. We assume a fully-radiative condition in this model 
(Ruffini et al. |2003[ l. The structures observed in the prompt 
emission of a GRB come from the inhomogeneities in this CBM, 
while in the standard fireball scenario ( Meszaros , 2006 ) they are 
caused by internal shocks. In this way we need few parameters 
for a complete description of a GRB: the dyadosphere energy 
E^gi, the baryon load B and the CBM density distribution, hcbm- 
In addition, we assume that there is spherical symmetry, and the 
energy released in the explosion £■,,„ is equal to the energy of 
the dyadosphere Ef^,. From this approach, to sum up, the GRB 
bolometric light curve will be composed of two main parts: the 
P-GRB and the extended afterglow. Their relative energetics and 
their observed time separation are functions of the parameters 
E'l^,, B, and ucbm- We want to stress that the emission of the P- 
GRB does not always coincide with what is called "prompt emis- 
sion" in the fireball scenario. Indeed, within the fireshell model, 
this prompt emission corresponds to the gamma-ray emission, 
which addresses not only the P-GRB, but also the peak of the 
extended afterglow. 

Instead of making use of the typical thermal spectrum, we in- 



-2.048: 45.057 s 



troduced a modified black body spectrum ( Patricelli et al. 2010 
|201 1) , given by 



dNy 
dVde 



(1) 



This way we can also reach an agreement with the most ener- 
getic GRBs {Eis„ > 10^^ erg). Furthermore, within the fireshell 
scenario we can naturally explain the hard-to-soft spectral varia- 
tion observed in the extended afterglow emission. As the Lorentz 
gamma factor F decreases with time, the observed effective tem- 
perature of the fireshell also decreases, making the peak of the 
emission take place at lower energies. This effect is amplified 
by the curvature effect of the EQTS ( [Bianco & Ruffinij |2005| l, 
which produces the observed time lag in the majority of the 
GRBs. 

We need to identify the P-GRB in the observed data so that 
we are able to determine the parameters El*, and B, via a trial and 
error procedure, and consequently the P-GRB energy Ep^cRB, 
the Lorentz gamma factor at the transparency y, the theoretically 
predicted temperature kT,ij, and the radius at the transparency 



(see Fig. 1 in (Ruffini et al. 2009 1). The observed temperature 



kTohs is related to the theoretically predicted temperature kT,i, 
through 



kT,h 
1+z' 



(2) 



5. Analysis of data and results 

To obtain the Fermi GBM light curve and spectrum in the band 
8-440 keV (see Fig.[2 1, we used the RMFIT program. We down- 
loaded the data from the gsfc websit^ We used the lightcurves 
corresponding to the second and fifth Nal detectors and the bO 
BGO detector We subtracted the background by fitting a cubic 
function from the intervals before and after the GRB (from 400 
s to 200 s before the GRB and from 180 s to 220 s after it). 



4>L02 + 





_02 + 
N'ALOS 




100 

Energy (kcV) 



ftp://legacy.gsfc.nasa.gov/fermi/data/gbm/bursts 



Fig. 4. Fit of the spectrum of episode 1, with a Band model (up- 
per panel) and a black body plus power-law model (lower panel). 



where we suppose there is no data. Then we proceeded with the 
time-resolved spectral analysis. 

To proceed with the fitting of the spectra, we defined first 
of all the time intervals we wanted to analyze: the first interval 
starts at the trigger time fo - and lasts 45 s, while the other 
starts at fo -H 45 s and lasts 44 s. For convenience, from now on 
we will refer to the first emission as episode 1 and the second 
emission as episode 2. For this source we considered two mod- 
els: the black body plus power-law model and the Band spectral 
model (Band et al.. 2003) . We first analyzed each of the events 
separately, as if they were two GRBs and then subdivided each 
of the two emissions in the light curve into two other parts: the 
one that we think would coiTespond to the P-GRB emission and 
the one that would coiTespond to the afterglow. The results from 
the spectral analysis are shown in Table [T] The fit of the spec- 
trum of the first episode with both models is shown in Fig. [4] 
while Fig. |5]shows the same fit for the second episode. 

6. Identification of the P-GRB 

6.1. Attempt for a single GRB scenario: the whoie emission 
as a single GRB 

The first step in our analysis was to attempt to interpret the whole 
emission as a single GRB, with episode 1 as the P-GRB. We 
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Table 1, Time-resolved specttal analysis of GRB 101023. 



Time [s] 


a 






^"'^""[keV] 


X- 


Norm 


kT[keW] 


r 


x' 


Norm '"' 


Norm 


0-44 


-1.3±0.8 


-1.94 


:0.2 


87±147 


0.98 


0.006±0.01 


14±6 


-1.7+0.1 


0.98 


0.0003+0.0004 


(4.1 +7.4) X 10-5 


45-89 


-0.9+0.1 


-2.02 


tO.l 


151±24 


1.09 


0.043+0.008 


26±1 


-1.58+0.03 


1.12 


0.0124+0.0006 


(4.2 +1.1) X 10"' 




Fig. 5. Fit of the spectrum of episode 2, with a Band model 
(upper panel) and a black body plus a power-law model (lower 
panel). Both models fit the entire energy range well, with a chi 
squared of 0.79 and 0.84, respectively. The data points have been 
grouped according a signal-to-noise ratio of N=10, and rebinned 
at higher energies in order to have better statistics and reduce the 
error bars. 



6.2. The identification of ttie P-GRB of the first episode 

Our second step in the analysis of this source was to attempt to 
interpret episodes 1 and 2 as two different GRBs. We first ana- 
lyzed episode 1 by taking two different possibilities into account; 

1 ) We considered a P-GRB that lasts 6 s and made the spec- 
tral analysis with XSPEC. We fitted a black body plus power-law 
model and found a black-body temperature of kT = 25.4 + 6.9 
keV with normalization factor normhhody - 0.9 + 0.5, a photon 
index of y = 2.2 ± 0.5 with nornipo - 30.9 + 35.3 and a re- 
duced chi squared of - 101. Considering that the P-GRB is 
the thermal component of the GRB, by using XSPEC we found a 
flux of 7.25 X 10"** erg/cm- /s in the range (8-5000) keV. Then we 
followed the same procedure for the whole of episode 1, fitting 
a cutoffpl model, and found a photon index of y = 1.16 ± 0.3, 
a cutoff energy of E cutoff = 73 + 27 keV, a normalization fac- 
tor of 2.9 + 2.4, a reduced chi squared value of )p_ - 1.08, and 
a flux of 1.626 x 10"^erg/cm^/s. Using formulaH we found a 
P-GRB energy of Ep^cRB = 9.56 x 10^*^ erg and a total energy of 
Ef*, = 1.625 X 10^2 erg, which gives a ratio Ep^crb = 5.9%£'f*,. 
With these values we performed the simulation with the numer- 
ical code and found a baryon load B = 8.5 x 10 "* and a pre- 
dicted temperature of kTii, - 128.82 keV, which is much higher 
than the one observed. Therefore, we concluded that the first 
6 s of emission cannot be the P-GRB of episode 1, at least in 
the fireshell scenario. 2) We considered the P-GRB under the 
threshold of the detector. We took the first 6 s before the trig- 
ger time as the P-GRB and supposed that it is well fitted by a 
Band model, with a flux of 10"^erg/cm^/s, which is comparable 
with the threshold of the detector. We derived a P-GRB energy 
of 10^" erg, which is the 0.9% of the total energy. For this ra- 
tio of the energies, we found with the numerical code a baryon 
load of Z? = 10 ~ and a predicted flux that is smaller than the 
detector threshold. This indicates that indeed this could be the 
P-GRB of the first emission, so that episode 1 could be a GRB, 
and we could be for the first time in the presence of a double 
GRB. However, in light of the results obtained from the analy- 
sis of GRB 090618 (Izzo et al.| |20ir) and taking into account 
that the value of the redshift has not been precisely determined, 
we decided to discard this result. Therefore, we conclude that 
episode 1 is not a GRB but another source whose origin is still 
unidentified. We come back to this interpretation later. 



performed a time-integrated analysis of the whole emission of 
episode 1, using a black body plus power-law model and a Band 
model. The results of this spectral analysis are shown in Table 
[T] We found a black-body temperature of kT - 14 + 6 keV with 
normalization factor normhhody = (4.1 + 7.4) x 10"^, a photon 
index of 7 = -!. 7 + 0.1 with normalization factor nornipo = (3 ± 
4) X 10 and ax^ ^ 0.98 for both spectral models. The P-GRB 
energy is Ep^qkb = 1.625 x lO*"^ erg and the isotropic energy 
Eiso - 4.03 X lO^^* erg, which gives a ratio Ep^qrbI Eiso — 0.04. 
This value in our simulations would imply a theoretically pre- 
dicted temperature of kT,h - 1 10.63 keV, which is by far much 
bigger than the observed one. Consequently, the first episode 
cannot be the P-GRB of the whole emission. 



6.3. Analysis of the second episode 

After the analysis of episode 1, we moved on to the analysis of 
episode 2. We followed the same steps taking the first 12 s of 
episode 2 as the possible P-GRB. We also fitted a black body 
plus power-law model to the whole P-GRB and found a black- 
body temperature of fer = 15.5+ 1.6 keV with normalization fac- 
tor normhhody = 1.26 ± 0.3, a photon index of 7 = 2.5 + 0.4 with 
normalization factor normpo = 141.79 and a^^ = 0.96. We com- 
puted a flux in flie band (260-5000) keV of 2.54x lO^^erg/cm^/s 
and a P-GRB energy of Ep^crb - 1.89 x 10^^ erg. By fitting a 
black body plus power-law model to the whole of episode 2 we 
found a flux in the band (8-5000) keV of 1.272 x lO^^erg/cm^/s 
and a total energy of E^^^ = 1.309 x lO^-' erg. The ratio is 
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Then we took the values of some parameters, the spec- 
trum, and response files from the XRT websit^ selected the 
part of interest and carried out an analysis making use of the 
program XSPEC. We fit the model wabs, which is the photo- 
electric absorption using Wisconsin cross sections (Morrison et 
T983] l: M(E) = exp{-nHa-{E)], where cr(£) is the photo- 



al 



electric cross section (not including Thomson scattering) and 
«H is the equivalent hydrogen column density, in units of 10^^ 
atoms/cm^. Once we knew these parameters, we fit the data 
with a power-law model, considering a phabs component re- 
lated to the intrinsic absorption. We obtained a value of n'^"' - 
0. 18±0.019x lO^^COT"^. Wkth this result, we put them in formula 
(1) of Grupe et al. ( 2007[ l paper: 



Time (s) 



log(l +z)<13- Q.5[log(l + ANh)], 

and we obtained an upper limit for the redshift of 3.8. 



(3) 



Fig. 6. Plot of the flux of the BB component vs time for the first 
12 s of episode 2. 



Ep-cRB - 0-9%£fo,- This same value is reached with the nu- 
merical code for a baryon load B = 7.6 x 10"^ and a predicted 
temperature of kT,/, - 14.02 keV, which after cosmological cor- 
rection gives 7.38 keV (assuming z=0.9, see next section), which 
is not in good agreement with the observed one, kTohs = 26 keV. 
Thus we conclude that the first 12 s of emission cannot be the P- 
GRB. 

To be more accurate, we performed the following procedure: 
as we know that the P-GRB consists of a black-body emission, 
we performed a detailed spectral analysis every Is with the Black 
body model to see the behavior of the black body component, i.e 
where the black body component dominates. That will indicate 
more precisely the time range and duration of the P-GRB. Table 
|2]shows the results of this analysis and Fig.|6]shows the behavior 
of the black body component with time. We see that in fact only 
the first 5 s of emission have a marked black body component, 
with a typical pulse shape. The emission that follows seems not 
to be related to the P-GRB, but to the afterglow. So we conclude 
that episode 2 is indeed a GRB and the first 5 s of emission are 
the P-GRB (see Section[8]l. 



7. Pseudo-redshift determination 

The redshift of this source is unknown, owing to the lack of data 
in the optical band. However, to constrain it, we employed three 
difl'erent methods, mentioned below. 



7.1. Method 1: nH column density 

We first tried to estimate the redshift making use of the method 
developed in Grupe et al. ( 2007 1 work, where the authors com- 
ment on the possible relation between the absorption column 
density in excess of the galactic absorption column density 
ANh = Nnjii - NH,gai and the redshift z. To do this, we con- 
sidered the galactic absorption component taken from [Kalberla 
|et al.| ( [2005| l, with the following values of the galactic coordi- 
nates of the GRB: I = 328.88, b = -38.88. We used the Lab 
Survey websit^and obtained the value of rin - 2.59 x 10^"cm"^ 
for the galactic H column density. 



7.2. Method 2: Amati relation 

We tried another method of constraining the redshift, making use 
of the Amati relation (Amati, 2006 ), shown in Fig.|7] This relates 
the isotropic energy Eiso emitted by a GRB to the peak energy 
in the rest frame Epj of its vFy electromagnetic spectrum (see 



Amati et al. 2009 and references therein). Eis,, is the isotropic- 



equivalent radiated energy, while Ep,, is the photon energy at 
which the time averaged vFy spectrum peaks. The analytical ex- 
pression of Eiso is 



47Tdf 



> boh 



(4) 



where is the luminosity distance, z is the redshift and S hoi is 
the bolometric fluence, related to the observed fluence in a given 
detection band (£■,„,„, -Emn.t) by 



S bol — S obs " 



Ji/(i+z) E4,{E)dE 



j^Z E<f>iE)dE 



(5) 



with (p the spectral model considered for the spectral data fit. The 
value of Epj is related to the peak energy Ej, in the observer's 
frame by 



Epj^Ep(l+z). 



(6) 



We started our analysis under the hypothesis that episode 2 
is a long GRB. We computed the values of Ep i and Eis,, for dif- 
ferent given values of z and plotted them in Fig. |7] We found 
that the Amati relation is fulfilled by episode 2 for 0.3 < z < 1.0. 
This interval has been calculated at Icr from the best fit from the 
Amati relation, in order to obtain a tighter interval around the 
best fit than with the previous method. 

7.3. Method 3: Empirical method for the pseudo-redshift 



We also tried an empirical method, following | Atteia| ( 20031) and 
|Pelangeon et al.|p006| l, which can be used as a redshift indicator. 
This method consists in determining a pseudo-redshift from the 
GRB spectral properties. Using the parameters from the Band 
model, namely the index of the low-energy power-law a and the 
break energy Eq, we can compute the value of the peak energy 



http://www.astro.uni-bonn.de/~webaiub/english/toolsJabsurvey.php 



http://www.swift.ac.uk/xrt_curves/ 
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Detailed spectral analysis of the P-GRB of episode 2, of 12s of duration, using a BB+po model and performed every 1 



r-OKK or episoue z 










Time Int 


Kl(KeV j 


norm lactoi 


rlux 0-44U KeV 


X 


051-052 


1.9±1.7 


0.9+2 


3.2525 xlO 


1.40 


03 2-05 J 


5±1 


1.3±0.3 




1.06 


053-054 


7±1 


1.2±0.3 


9.9689x10 


0.99 


054-055 


10±2 


1.2±0.3 


9.8285x10"** 


1.17 


055-056 


7±1 


1.6±0.3 


1.3217x10-' 


0.96 


056-057 


10±1 


2.1±0.4 


1.7721x10-' 


1.42 


057-058 


10±1 


1.7±0.4 


1.4245x10"' 


0.96 


058-059 


11±1 


2.1±0.4 


1.7738x10-' 


1.16 


059-060 


10±1 


2.6±0.4 


2.1844x10-' 


1.38 


060-061 


10±1 


1.8±0.3 


1.4976x10-' 


1.51 


061-062 


9±1 


1.8±0.3 


1.5193x10-' 


1.18 


062-063 


14±2 


1.6±0.4 


1.3462x10-' 


1.74 



En peak(keV) 



1000 



100 - 




10 - 



En iso(erg) 



Fig. 7. Plot of the relation between Epj and Ejso for the second 
episode of GRB 101023, considering different values of the red- 
shift. It can be seen that the plot lies within Icr for the range z- 
0.3 -z= 1.0. 




Jims since burst trigger (s) 

Fig. 8. Fit of the second major pulse of the light curve of GRB 
101023. 



of the vFy spectrum, as Ep - Eq{2 + a). Then, we define the 
isotropic-equivalent number of photons in a GRB, Ny, as the 
number of photons below the break, integrated from 00 to 
Ep/2. If we also know the Tgo, we define the redshift indicator 



X 



/Tgo 



(7) 



From a sample of 17 GRBs with known redshift reported in 
Atteia| (|2003)l we compute the theoretical evolution of X with the 
redshift z, that is X = f(z). Then we invert the relation to derive 
a pseudo-redshift from the value of X. That way we obtain the 
pseudo-redshift as z = /-'(X), for the GRB of interest. 

We applied this treatment to episode 2 of GRB 101023, 
introducing the spectral parameters from the Band model on 
the Cosmos website0 and obtained a value for the redshift of 
z - 0.9 + 0.084. It is important to mention here that this error is a 



statistical one, while the systematic error is much bigger (jAtteia 
[2003] IPelangeon et aLl |20D6l |PeIangeon| [2008| , of a factorof^ 
1.5, i.e., z = 0.9!^ 5^ 

This result agrees with the redshift range found from the 
Amati relation for episode 2 and is also consistent with the upper 
limit determined with method 1 . 

http://cosmos.ast.obs-mip.fr/projet/v2/fast_computation.html 



8. Simulation of the light curve and spectrum 

To simulate the light curve we made use of a numerical code 
called GRBsim. This numerical code simulates a GRB emis- 
sion by solving the fireshell equations of motion, taking the ef- 



fect of the EQuiTemporal Surfaces (EQTS, .Bianco & Ruffini 



2005) into account. We made the simulation for episode 2. We 
found, at the transparency point, a value of the laboratory ra- 
dius of 1.34 X 10'"* cm, a theoretically predicted temperature 
that after cosmological correction gives kT,i, - 13.26 keV, a 
Lorentz Gamma factor of F = 260.48, a P-GRB laboratory en- 
ergy of 2.51 X 10^' erg and a P-GRB observed temperature of 
28.43keV. We adopted a value for the dyadosphere energy of 
Ef* = 1.8 X 10" erg and a baryon loading of B = 3.8 x 10-^ 
The simulated light curve and spectrum of episode 2 are shown 
in Figs.[8]and[9j respectively. 

Figure [To] shows the fitted spectrum with different models. 
We took the data points from the Nal n2 and BGO bO detectors 
together We note there is a good agreement between both fits, 
in the low and medium energy range. At high energies, the spec- 
trum follows a power-law behavior, which cannot be reproduced 
by the modified black body model due to the exponential cutoff. 
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Spectrum Fermi GBM (8-440 keV) i 
Simulation ot the spectrum 



ii. I 




Energy(keV) 



Fig. 9. Fit of the spectrum of episode 2. 




100 200 



Eiierg\'(keV) 



Fig. 10. Fit of the spectrum of episode 2. The green dotted lines 
represent the fit of a black body and a power-law components, 
separately. The red line is the sum of them, calculated with 
XSPEC (BB-i-po). The blue line is the fit with the modified black 
body spectrum given in Eq.(|T]l, calculated with the GRBsim nu- 
merical code. 

9. Analysis of the first episode 

To analyze episode 1 more into detail, in order to identify the 
nature of this phenomenon, we plotted the temperature of the 
black body component as a function of time, for the first 20 s of 
emission (see Fig.[TT|. We note a strong evolution in the first 20 s 
of emission which, according tolRyde ( 2004) 1 can be reproduced 
by a broken power-law behavior, with a = -0.47 + 0.34 and 
yS - -1.48± 1.13 being the indices of the first and second power 
law, respectively. We also plotted the radius of the most external 
shell with time (see Fig. 
radius can be written as 



12 1. Following Izzo et al. (2011 1, the 



RDF 



(l+z)2' 



(8) 



where R 



(pobs/i^^crT^^i^^) is a parameter, D the luminosity 
distance, F the Lorentz factor, and <poi,s the observed flux. We can 
see that the radius remains almost constant (in fact it increases, 
but only slightly). From this it is possible to see that the plasma 



k T(keV) 
5fl 



30 



+ 



ft 



1.5 2.0 



3.0 



5.0 7.0 



10.0 



15.0 20.0 



Fig. 11. Evolution of the observed temperature kT of the BE 
component. The blue line corresponds to a broken power-law 
fit. The indices of the first and second power laws are ff = 
-0.47 + 0.34 and yS = -1.48 ± 1.13, respectively. The break 
occurs at 1 1 s after the trigger time. 

Radius(km) 
IxlO" 



5x10' 

2x10^ 
1x10' 
5x10* 

2x10"* 



+ 4+ 



t. 



Time{s) 

1.5 2.0 3.0 5.0 7.0 10.0 15.0 20.0 

Fig. 12. Evolution of the radius of the first episode progenitor. 



is expanding at nonrelativistic velocities. According to the work 
( 2011| ), there is an expansion phase of the 
e the iron core suffers a contraction. This is 



of lArnett & Meakin 



boundary layers, whili 
due to the presence of strong waves originated while the different 
shells of the progenitor mix during the collapse phase. This fact 
confirms the non-GRB nature for the first episode. 

9.1. The X-ray afterglow as a possible redshift estimate ? 

We have seen that GRB 090618 and GRB 101023 share similar 
properties. They seem to be composed of two different emission 
episodes, the first being connected to a quasi-thermal process 
before the collapse of the core, while the second is the canonical 
GRB (see |Ruffini et al.[ 2010a d). 

Anyway, if both GRBs were created originated by the same 
physical mechanism and since the energetics are very similar, 
considering the value z = 0.9 for GRB 101023, we can expect 
similar luminosity behavior for the X-ray afterglow. Although 
we have not yet developed a theory for this late afterglow emis- 
sion, we attempted a simple test that compared the observed X- 
ray afterglow of both GRBs as if they were located at the same 
redshift. Since there are different spectral components in the 
GRB X-ray afterglow, we built the pseudo-redshift light curves 
for both these different emissions. Thanks to the Swift-XRT ob- 
servations, we know that the early X-ray afterglow of both GRBs 
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Fig. 13. The fit of the X-ray afterglow of GRB 090618 (up- 
per panel) and GRB 101023 (lower panel) with the model of 
WilUngale et aL]([2007l). 



shows a canonical behavior, where the emission can be divided 
in three distinct parts (Nousek 2006): 1) a first very steep de- 
cay, associated with the late prompt emission; 2) a shallower de- 
cay, the plateau; 3) a final steeper decay. At first, we determined 
for GRB 090618 and GRB 101023 these three time intervals by 
using the phenomenological function introduced in the work of 
IWilUngale et aL| ( |2007l ): 



10' 




Fig. 14. The X-ray afterglow of GRB 090618 (blue data) as if it 
was observed at redshift z - 0.9 (see text). The X-ray afterglow 
of GRB 101023 is also shown as comparison (red data). Data on 
GRB 101023 are missing between ~ 200 s and 3550 s. Where 
data are present, the superposition is striking. 

assumption that the spectrum of each time interval is best fit- 
ted by a simple power-law model. This assumption is supported 
by the hypothesis that the X-ray afterglow comes from a syn- 



chrotron emission mechanism Sari et al. ( 1999 1, whose spectral 
emission is represented by a simple power law function. Then, 
we extrapolated the emission of the afterglow of GRB 090618 
in the (0.57 - 19) keV energy range by considering the ratio be- 
tween the number of photon counts in both energy ranges. This 
value corresponds to a conversion factor, which we consider for 
scaling the intensity of the light curve. We finally amplified, by 
a term (1 + Zi()io23)/(l + -Zoweis)- the time interval of emission 
of GRB 090618, obtaining as a final result the afterglow light 
curve of GRB 090618 as if it was observed by XRT at redshift 
0.9, see Fig. 14 It is, most remarkably, a perfect superposition of 



f{t)^{Fcexp\a, 



exp 



(^),r<r.;F.(l-)"^x,(^) 



which represents the transition from an exponential regime to a 
power law. This transition occurs at the point (Tc, Fc) where the 
two functional sections have the same value and gradient. The 
parameter determines both the time constant of the exponential 
decay and the temporal decay index of the power law, while the 
tc parameter marks the initial rise. The maximum flux occurs at 
t - (tcTc/acy^^. We fit the afterglow data of the two GRBs with 
this model, and the results of our fits are shown in Fig.[T3] 

After the determination of these three time intervals, we built 
the X-ray light curve of GRB 090618 as if it was observed at 
redshift z - 0.9, which is our estimate for the redshift of GRB 
101023. The Swift-XRT (which operates in the (0.3 - 10) keV 



the light curve emission of both GRBs. This evidence delineates 
three important aspects: 

- the X-ray afterglow of both GRBs clearly confirms a com- 
mon physical mechanism for these GRBs; 

- there is ample convergence and redundancy with different 
methods of determining a value of redshift z = 0.9 for GRB 
101023. There has also been the unexpected result pointing 
to the late afterglow as a possibly independent redshift esti- 
filia!0^;(9) 

the redshift of GRB 101023 derived by the superposition of 
the two afterglow curves is consistent with the value of z = 
0.9, which we have found before. 



This last point led us to do another analysis consisting in 
the redshift-translation of the X-ray afterglow of GRB 090618 
considering different values for the redshift. Following the same 
procedure and considering five different values for the redshift, 
z = (0.4,0.6,1.2,2,3), we see that the X-ray emission of GRB 
101023 is compatible with the X-ray afterglow of GRB 090618 
as if it bursted between z = 0.6 and z - 1.2, see Fig. 15 Then 



energy range) light curve of GRB 090618 (Evans et al. 2007 



2009[ l corresponds to the emission in the rest frame at z = 0.54 in 



the energy range (0.462 - 15.4) keV, while for GRB 101023 the 
XRT window corresponds to the range (0.57 - 19) keV. We must 
obtain the emission of GRB 090618 in this last energy range, 
in order to compare the two light curves. At first we made the 



we conclude that our estimate for the redshift of GRB 101023 of 
z = 0.9 is very reliable. 

10. Conclusions 

GRB 101023 is a very interesting source for the following rea- 
sons. 

1) We find a striking similarity between GRB 101023 and 
GRB 090618, as can be seen from the light curves. Following 
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io-'2 



BOO 1000 lo-" 10- 




Fig.15. The X-ray afterglow of GRB 090618 as if it was ob- 
served at different redshifts z - (0.4,0.6, 1.2,2,3), where each 
color corresponds to a different redshift. The X-ray afterglow of 
GRB 101023 is also shown for comparison (red data). 



the study of GRB 090618, we divided the emission into two 
episodes: episode 1, which lasts 45 s, presents a smooth emission 
without spikes that decays slowly with time. Episode 2, of 44 s 
of duration, presents a spiky structure, composed of a short and 
faint peak at the beginning, followed by several intense bumps, 
after which there is a fast decay with time. Episode 2 has all the 
characteristics of a canonical long GRB. 

2) We performed a time-resolved analysis of episode 1 . We 
fitted a black body plus a power-law model and plotted the evolu- 
tion of the black body component with time. The observed tem- 
perature decreases during the first 20 s following a broken power 
law: the first with index a - -0.47 + 0.34 and the second with 
index ;6 = -1.48 + 1.13, see Sec. 9. This behavior is very similar 
to GRB 090618. 

3) In the absence of a direct measurement of the redshift to 
the source, we have inferred it from several empirical methods. 
First, following the work of Grupe et al. ( 2007 1, which consid- 
ers the hydrogen equivalent column density in the direction of 
the source, we obtained an upper limit of z < 3.8. Then we per- 
formed a spectral analysis to episode 2, fitting a Band model. 
From the peak energy Epeak and using the Amati relation un- 
der the hypothesis that episode 2 is a canonical long GRB, we 
constrain the value of the redshift to be between 0.3 and 1.0. 
Finally, using the parameters of the Band model and following 
the work of |Atteia| ( |2003[ ), we determine a value of the redshift 
of z = 0.9 + O.QMistat.) + 0.2{sys.). The thi-ee methods are con- 
sistent, so we assumed for the redshift of this source z - 0.9. 

4) From the knowledge of the redshift of the source, we have 
analyzed episode 2 within the fireshell model. We determined a 
total energy Eiso - 1 .79 x 10^^ erg and a P-GRB energy of 2.5 1 x 
10^' erg, which we used to simulate the light curve and spectrum 
with the numerical code GRBsim. We find a baryon load B - 
3.8x 10"^ and, at the transparency point, a value of the laboratory 
radius of 1.34 x 10'^ cm, a theoretically predicted temperature of 
kT,ij - 13.26 keV (after cosmological correction) and a Lorentz 
gamma factor of F = 260.48, confirming that episode 2 is indeed 
a canonical GRB. 

5) From the knowledge of the redshift, we can also evalu- 
ate the flux emitted by episode 1, and from the observed black 
body temperature, infer the radius of the black body emitter and 
its variation with time, see Fig. 12. We saw that it increases 
during the first 20 s of emission, with a velocity ~ 1.5 x 10"* 
km/s. In analogy with GRB 090618, we concluded that episode 
1 originates in the last phases of gravitational collapse of a stel- 



lar core, just prior to the collapse to a black hole. We call this 
core a "proto-black hole" (! Ruffini et al.[ 2010a). Immediately 
afterwards, the collapse occurs and the GRB is emitted (episode 
2). 

6) Finally, we performed the following test. Owing to the 
similai-ities between GRB 101023 and GRB 090618 regarding 
morphology and energetics, we expect them to be created by the 
same physical mechanism, so we compared the late observed 
X-ray afterglow of both GRBs as if they were located at the 
same redshift; i.e, we built the light curve of GRB 090618 (of 
z = 0.54) as if it had redshift z - 0.9, exti-apolating it to the XRT 
energy window of GRB 101023. We found a surprising perfect 
superposition of the light curves for z=0.9, receiving a further 
confirmation of the correctness of the cosmological redshift de- 
termination. The same procedure for the redshift determination 
will be repeated for sources with a spectroscopical-determined 
redshift, as a further check of our proposal. This result points to 
a possible use of the late afterglow as a distance indicator. 

We concluded that GRB 101023 and GRB 090618 have 
striking analogies and are members of a specific new family of 
GRBs developing out of a single core collapse. It is also appro- 
priate to remark that this new kind of source does not present any 
GeV emission. The existence of precise scaling laws between 
these two sources opens a new window on the use of GRBs as 
distance indicators. We will go on to identify additional sources 
belonging to this family. This new paradigm is also being ap- 
plied to sources at very high redshift to see how the absence of 
a signal under the threshold can affect the theoretical interpreta- 
tion. We are also considering the possibility that proto-neutron 
stars in addition to proto-black holes may exist in the case of su- 
pemovae or hypemovae. Particularly interesting in this respect 
is the work of So derberg et al.| (j2008) showing the X-ray emis- 
sion prior to SN events, which may relate the observed X-ray 
emission prior to SN 2008D to episode 1 in GRB 090618 and 
GRB 101023. In this sense we are revisiting our considerations 



of GRB 980425 (see e. g.|Fraschet ti et al. 2 004, Fraschetti etal] 
2005 ; Ruffini et al. 200 4a||2007[|B ernardin i et ai.||2008| l, as""well 
as of G RB 030329 (|Bernard ini et al. 2004', '20()5b') and GRB 
031203 ( |Bernardini et al.| |2005a, Ruflini et al. 2007 2008). 
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